Critical scaling in random-field systems: 2 or 3 independent exponents? 
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We show that the critical scaling behavior of random-field systems with short-range interactions 
and disorder correlations cannot be described in general by only two independent exponents, con- 
trary to previous claims. This conclusion is based on a theoretical description of the whole {d, N) 
domain of the rf-dimensional random-field 0{N) model and points to the role of rare events that 
are overlooked by the proposed derivations of two-exponent scaling. Quite strikingly, however, the 
numerical estimates of the critical exponents of the random field Ising model are extremely close 
to the predictions of the two-exponent scaling, so that the issue cannot be decided on the basis of 
numerical simulations. 

PACS numbers: ll.lO.Hi, 75.40.Cx 



I. INTRODUCTION 

The critical behavior of models in the presence of a 
quenched random field has attracted a lot of attention 
over the past decades.-^i^ Beyond the experimental inter- 
est, such models provide a rich playground to investigate 
the infiuence of quenched disorder on the long-distance 
properties of a system. Their equilibrium critical point, 
separating a disordered, paramagnetic phase from an or- 
dered, ferromagnetic one, is known to be controlled by 
a zero-temperature fixed point, with temperature being 
an irrelevant variable in the renormalization-group sense. 
Among the new features brought by such a fixed point 
is the violation of the "hyperscaling relation" between 
critical exponents. The space dimension d that appears 
in this relation must be corrected by the exponent 9 
describing the flow of the renormalized temperature to 
zero, thereby leading to an apparent dimensional reduc- 
tion from d to d — Oik'^ 

Whereas phenomenological theories take 9 as an in- 
dependent exponent ^i^ which implies that equilibrium 
scaling behavior is described by three independent ex- 
ponents in place of the usual two-exponent scaling for 
finite-temperature fixed points, two different approaches 
have claimed that 9 is actually fixed. One is the 
dimensional-reduction prediction, based either on pertur- 
bation theory^^^ or on the Parisi-Sourlas supersymmetric 
formulation J. It states that 9 — 2 and that the exponents 
of the random-field model are those of the correspond- 
ing pure model in two dimension less. This prediction 
has been rigorously proven to be wrong in low enough 
dimension,^' ^"^ and we have recently provided a complete 
theoretical explanation of dimensional reduction and its 
breakdown through a nonperturbative functional renor- 
malization groupJ^"— 

Another line of argument has been put forward by 



Schwartz and coworkers. ^^'^^ The claim is that 9 = 2 — r], 
with r] the anomalous dimension of the (order parameter) 
field, so that scaling is described by only two independent 
exponents, e.g., rj and the correlation length exponent v. 
The derivation involves general manipulations and the 
result is supposed to hold for the Ising as well as the 
continuous version with 0{N) symmetry. It is also sup- 
ported by heuristic considerations J^J^i 

The problem with the above derivations is that they 
rely on formal manipulations that are blind to the pres- 
ence of rare events or rare regions, such as avalanches 
or droplets, or that overlook the presence of multiple 
metastable states. This is known to be the reason for the 
failure of the simple supersymmetric formulation . ^'^'^^i^^ 
This also casts some doubt on the general validity of the 
order-of-magnitude estimates of the relative strength of 
the fiuctuations carried out in Schwartz's arguments J^ii^ 

How to conclude on the validity of the two-exponent 
scenario proposed by Schwartz and coworkers? The ques- 
tion is actually quite subtle as the description is exact 
at the lower and the higher critical dimensions of the 
random-field models and it appears to be numerically 
very well verified in computer simulations (including a 
recent extensive one^^), high-temperature expansions or 
approximate renormalization-group treatments, mostly 
in d = 3. Our contention is that it is impossible to answer 
the question on a pure numerical ground as the error bars 
will always blur the conclusion. On the contrary, with the 
help of the functional renormalization group (FRG), the 
problem can be studied for continuous values of the di- 
mension d and of the number of components N of the 
order-parameter field. We can then show with no ambi- 
guity that the two-exponent scenario cannot be right in 
general, which then reduces the corresponding prediction 
for the random-field Ising model in d = 3 to a very, or 
even extremely, '^^ good but by no means exact, result. 



We have already made this point in earher worki^ii^ but 
in this paper, we sharpen the arguments and present ex- 
tended calculations. 

The rest of the paper is organized as follows. In sec- 
tion II, we present the the random field 0{N) model 
[RFO(A^)M] and the scaling behavior around its criti- 
cal point and in section III we summarize and discuss 
the predictions of the two-exponent scaling scenario. In 
section IV, we recall the results of the FRG for the 
RFO(iV)M with A^ > 1 (continuous symmetry) near the 
lower critical dimension for ferromagnetism, in d = 4-|-e. 
We show that at order e and e^, the relation proposed 
by Schwartz and coworkers is unambiguously violated. 
In the next section, we consider the random field Ising 
model (RFIM) via a nonperturbative version of the FRG 
that we have previously developed. Here too, we show 
that the two-exponent scenario cannot be right for all di- 
mensions between 2 and 6. We focus in particular on the 
dimensions around the critical dimension cIdr — 5.1 that 
marks the lower limit of existence of the supersymmetric 
fixed point associated with the d — > d — 2 dimensional 
reduction. Finally, we give some concluding remarks in 
section V. 



II. THE CRITICAL BEHAVIOR OF RANDOM 
FIELD SYSTEMS 



The long-distance behavior of the random field 0{N) 
model is described by the following Hamiltonian or bare 
action, 



s[vM 



-\dcp{x)f + UB{\<pix)\') ~Hx)M^) 



(1) 

where J = J d'^x, (p{x) is an N component field, 
UBi\<f\^) = (r/2)|¥>|2 + (u/4!)(|¥'|2)2, and h{x) is a ran- 
dom source (a random magnetic field in the language of 
magnetic systems) with zero mean and a variance 



ht^{x)h-{y)=ABS^j'^''\x-y) 



(2) 



where //, z^ = 1, • • • , N and an overline denotes an average 
over the random field. An ultra-violet (UV) momentum 
cutoff A, associated with an inverse microscopic length 
scale such as a lattice spacing, is also implicitly consid- 
ered. We focus here on the short-range version defined 
above. We shall briefiy comment on the long-range ver- 
sion in the conclusion. 

Due to the presence of the random field, one needs 
to consider two different types of pair correlation func- 
tions of the y ) field: the so-called connec ted one, 

GHI^ni^ - V) 

disconnected 



= {ipt^{x)Lp''{y)) - {ipf^{x)){ip''{y)) and the 

^ne, G^f,(x - y) = {ip^^ (x)) {ip- (y)) - 

{ipi^{x)) {ip'^(y)). At the critical point Tc, the two cor- 



relation functions behave as 
Gconi^ ^ y) ^ 

GZi^-y)- 



TSn,, 



\x-y\'^-^+f^ ' 



(3) 



where rj is the usual anomalous dimension of the field 
and ^ is a priori a new exponent. Accordingly, one can 
also define two types of susceptibilities that diverge as 
one approaches the critical point from above as 



Xcon — I Gcon\X) ^ (T — Tcj 

Xdts = / Gdis{x) ^ {T - T^y , 



(4) 



the former one being the usual magnetic susceptibility. 
[The component indices have been dropped as the func- 
tions are then proportional to 6^^, as in Eq. ([3|).] The 
exponents 7 and 7 are related to rj and fj via 7 = (2~ri)u 
and 7 — (4 — 77)1/, with v the correlation length exponent. 
The renormalized temperature is irrelevant at the fixed 
point controlling the critical behaviori^ and it flows to 
zero with an exponent 9. As a result, the hyperscaling 
relation has an unusual form. 



{d-B)v , 



(5) 



where a is the specific-heat exponent. The exponents 0, 
ry, and fj are related through 



= 2 



+ V-V , 



(6) 



so that the scaling around the critical point is a priori 
described by three independent exponents, e.g., 77, 1/, and 
6 or fj. 



III. THE TWO-EXPONENT SCALING 
DESCRIPTION 

The d — > d — 2 dimensional reduction predicts that 
9 = 2, i.e. fj — rj, and furthermore that all the critical 
exponents are given by those of the pure model in (d — 
2). On the other hand, the two-exponent scenario put 
forward by Schwartz and coworker a^^'^^ states that the 
exponents obey the following relations: 



= 2 — rj , rj = 2ri . 



(7) 



The derivation actually also implies that the discon- 
nected and connected correlation functions in Fourier 
space are related through 



Gdis{q) = -j^ Gcon{qY 



(8) 



at criticality and when g — > 0. Eq. (JH) implies that the 
second cumulant of the random field is not renormalized 
and stays fixed to its bare value As . 



A stronger claim was initially made by Schwartz/^ who 
suggested that all the exponents of a random-field sys- 
tem in dimension d are the same as those of its pure 
counterpart in a reduced dimension d — 2 + ri[d). This 
prediction was however soon shown by Bray and Moore^^ 
to be already wrong for the exponent i' in the RFIM near 
its lower critical dimension 2 at first order in e = d — 2. 
To the best of our knowledge, it was subsequently aban- 
doned. 

It is also worth mentioning that various heuristic argu- 
ments have been used to derive the above two-exponent 
scaling»i^iii For instance, one such argument states that 
the magnetization per spin in a finite-size system of linear 
size L (we consider here the RFIM for simplicity) , which 
scales as L~('^~4+^)/2^ jg given for a typical realization of 
the random field by the average magnetic susceptibility 
Xcon,L ^ i^^'' times the mean random field h^ which 
scales as L^'^'^: 



rriL ^ L 



-(d-i+fi)/2 



Xcon,L hi 



L^-VL-d/2 ^ 



(9) 



It immediately results that fj = 2r]. 

In the case of the RFIM, the two-exponent scenario 
is exact at and near the lower critical dimension in first 
order in e = d — 2, as shown in Ref. J2C|, and it is 
also somewhat trivially expected at the upper dimension 
d = 6 and in first order in e = 6 — d as both rj and f] 
are zero and the mean-field result 9 — 2 still applies. In 
between, numerical studies via high-temperature expan- 
sions in d = 3, 4, 5^ and computer simulations in d = 3, 
confirm that 7 is very close to 27, and 77 to 2?], certainly 
within the accuracy of the methods. 

However, as mentioned in the Introduction, a funda- 
mental problem remains that the proposed derivations of 
the two-exponent scaling involve, beyond formal manip- 
ulations, estimates of the relative order of magnitude of 
the fluctuations that essentially rely on factorization ap- 
proximations and the central-limit theorem in the limit 
of large system size (see e.g. the above heuristic argu- 
ment). All of the derivations are therefore blind to rare 
events, rare regions or rare samples, which, precisely, 
have been shown to be crucial in disordered systems near 
zero-temperature fixed points 4i^ii2r2ii22ri2^. 

In the absence of rigorous derivations, numerical ev- 
idence, no matter how good, is insufficient to establish 
the validity of the scenario, because of unavoidable un- 
certainty. Quite the contrary, we show below that the 
two-exponent description, which is claimed to apply to 
random-field systems below their upper critical dimen- 
sion irrespective of the number of components N and of 
the dimension d, cannot be valid in general and, as a 
consequence, has no rigorous foundations at this point. 



IV. THE RFO(N)M NEAR d = 4 

Near the lower critical dimension for ferromagnetism, 
which is equal to d = 4 in this case, the long-distance 
physics of the RFO(iV)M is captured by a nonlinear a 



model that in turn can be studied through a perturbative 
but functional RG^ The resulting FRG flow equations 
in d = 4+e have been studied at one i^i^^i^^ g^^^j two ^^1^^ 
loop order. 

The dimensionless second cumulant of the renormal- 
ized random field is a function A(z) (where z is the co- 
sine of the angle between fields in two different copies of 
the system^^'^^) and it obeys the following RG equation 
at one loop: 



kdkAk{z) = eAfc(z) 



zAfc(z)2 + (iV-3)Afe(l)Afc(z) 



+ (A^ - 3 + 4z2)Afc(z)Ai,(z) -iN + l)zAkil)Kiz) 
- z(l - z^)A,{z)AUz) + (1 - z^)A,{l)A'az) 

~ 3z(l - z^)A',izr + (1 - zYKiz)AUz) , 

(10) 
where k is the running infrared (IR) momentum cutoff, 
a prime denotes a derivative, and Afc(z) [which, up to 
a multiplicative constant, was denoted R'{z) in previous 
publications] is of order e near the fixed point. One can 
moreover define two running exponents rjk and f]k as fol- 
lows: 



77fc - Afe(l), fik = -e+{N- l)Afe(l) 



(11) 



They converge to the fixed point values 77, fj when fc — >■ 0. 
The cor resp onding equations at two-loop order are given 
in Ref. |23l and are not reproduced here. 

A numerical and analytical investigation of these FRG 
equationsiSiS^iS^ shows that above a critical value of the 
number of components, Ndr = 18 — ^e, there exists a 
fixed point that corresponds to the d — > d— 2 dimensional 
reduction, with f\ = -q = jf^ — (jv~ 2)^ ■ ^^ove a slightly 
higher value, N^usp = 2(4 + 3%/3) - 3(2 -I- 3V3)e/2,^ this 
fixed point is stable and describes the critical behavior 
of the RFO(A^)M (see also Ref. [39]). In this domain of 
N , critical scaling is therefore described hy fj — rj and 
9 = 2, which contradicts the predictions of Schwartz and 
coworkers. -i^ii^ 

Below Ncusp, the stable fixed point is now character- 
ized by a nonanalyticity in the functional dependence of 
the renormalized disorder cumulant [A^(z) — A^(l) ^ 
Vl ~ z when 2; — ^ 1] that is strong enough to break the 
dimensional-reduction prediction, with therefore fj > rj. 
However, the possibility that this nonanalytic "cuspy" 
fixed point is characterized by the equality 77 = 277 is in- 
validated by the results. We display the ratio (2ry — fj)/ri 
as a function of A^ at one-loop order in Fig. 1. (The re- 
sults are confirmed at two-loop order.) The ratio is equal 
to 1 above Ncusp — 18.3923 and decreases continuously 
as A^ decreases. At the endpoint value at which A*(l) 
diverges, Nc = 2.8347- • • , it reaches a strictly positive 
value of 3 — A^c — 0.1653 • • • . We stress that the out- 
put can be obtained with an arbitrary precision, which 
is quite different from simulations. 




2.8347 



FIG. 1: Ratio (2r; — ff)/r] as a function of A'' for the 
RFO(A'')M in d = 4 + e at one-loop order. Above Ncusp = 
2(4 + 3\/3) — 18.3923 • • ■ the stable fixed point is associated 
with the d — >■ d — 2 dimensional reduction and {2ri — fj)/rj = 1. 
Below Ndr ~ 18 (marked by the leftmost vertical dashed 
line), the critical behavior is described by a "cuspy" fixed 
point and dimensional reduction breaks down but (277 — f])/ri 
is generically > 0. In a small interval Ndr. < N < Ncusp the 
two fixed points coexist, but the dimensional-reduction one is 
unstable to a nonanalytic, "cuspy" perturbation. 



V. THE RFIM NEAR cIdr i; 5.1 

The main result proving that the two-exponent scaling 
is generically inexact in the short-range RFIM is that 
there exists a range of dimension below the upper critical 
one, d = 6, for which the d — )■ d— 2 dimensional reduction 
is valid and therefore for which fj = rj ^ 2r] and 9 = 
2 y^ 2 — T]. The arguments establishing this result are as 
follows: 

(i) The fixed point associated with dimensional reduc- 
tion proceeds continuously from the Gaussian fixed point 
in d = 6 as one lowers d. It is characterized by the ab- 
sence of a linear cusp in the functional dependence of the 
cumulants of the renormalized random field. (It can in- 
deed be shown by the FRG that only such a cusp can 
lead to a breakdown of dimensional reductioUi ^^'^'^ ) 

(ii) That this dimensional-reduction fixed point is sta- 
ble for some interval of d below 6 can furthermore be seen 
by studying, in addition to all possible "cuspless" per- 
turbations that indeed prove to be irrelevant, the eigen- 
value associated with a nonanalytic, "cuspy" perturba- 
tion around this fixed point. This can be done exactly 
in d = 6 and can be numerically obtained below 6. The 
perturbation is then found to be irrelevant at and in the 
vicinity of d ~ 6: the eigenvalue is equal to 1 in d = 6 
and continuously decreases as one lower d slightly below 

(iii) By means of a nonperturbative truncation of the 
exact hierarchy of FRG equations (NP-FRG) for the cu- 
mulants of the renormalized disorder, we have located 
the limit of existence of the dimensional-reduction fixed 
point at doR ~ 5.1ji2, There is thus a nonzero range of di- 
mension above d_D_R and below 6 in which fy = 77 7^ and 



9 = 2, which contradicts the claim of the two-exponent 
scaling scenario. 

In addition, we have investigated in more detail the 
"cuspy" fixed point corresponding to a breakdown of the 
d — > d — 2 dimensional reduction below d^j^. Our pur- 
pose was to show that this fixed point is described by 
three independent exponents in general. The NP-FRG 
equations that must be numerically solved are given in 
Ref. [1J|, together with technical comments, and this is 
not reproduced here. 

We focus on the vicinity of doB., which is where the vi- 
olation of the proposed relation 77 = 277 is unambiguous. 
We display in Fig. 2 for illustration the RG flow of the 
running exponents 77^ and 77^ (for a crisper visualization, 
we actually plot the combination 27^^ — f/fc) in d = 5 as a 
function of log(fc/A), where A is the UV cutoff and k the 
running IR scale. When starting from initial conditions 
at the UV scale that are analytic (as the bare action), rjk 
and fjk first equal each other until one reaches a scale, the 
so-called "Larkin scale" ji^i^iiSS ^^ which a strong enough 
nonanalyticity appears in the functional dependence of 
the second cumulant of the renormalized random field 
and the two exponents 77^ and fjk start to deviate. How- 
ever, as can be seen, d = 5 is close to doR, so that the 
difference between the exponents remains small at the 
fixed point (fc — >■ 0) and 27; — ^ is strictly larger than 0. 

We also plot in Fig. 3 the fixed-point value of the 
ratio (277 — 77)777 as a function of d for the RFIM in the 
vicinity of djjji from the NP-FRG equations:^ The crux 
of the present demonstration is not the actual values of 
the exponents (which, however, are always within 10% 
of the best estimates), but the fact that, around doR, tj 
and 77 are equal or close to each other so that (277 — ??)/?7 
stays near 1 and unambiguously violates the prediction 
of the two-exponent scenario. This is clearly seen in Fig. 
3. 

This conclusion is similar to that reached above for 
the RFO(Af)M with continuous symmetry near d = 4 
and it may actually be extended through the NP-FRG 
to the whole (N, d) domain^- (note that the approximate 
nonperturbative truncation of the exact FRG equations 
gives back the exact perturbative FRG treatment in d = 
4 + e at one loop) . 



VI. CONCLUSION 

In this paper, we have challenged the two-exponent 
scaling scenario of the critical behavior of random-field 
systems proposed by Schwartz and coworkerSi^^'^^ Due 
to the generality of its derivations, this scenario is sup- 
posed to apply to all random field models, whether in 
the Ising or in the 0{N) version, and in all dimensions 
d between the lower and the upper critical ones. We 
have however clearly proven that the predictions, most 
notably the relation between exponents 77 = 277, cannot 
be exact for generic random-field systemSf^ Our proof 
is beyond error bars as it involves the fact that there is 
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FIG. 2: NP-FRG flow of the running exponent i]k and of 
2?7fe —fjk as a function of f = log(fc/A) for the RFIM in d = 5. 
The arrows indicate the flow toward the infrared {k — >■ 0). 
Note that one reaches i ~ — 19, which means a length scale 
bigger than 10* times the microscopic one. One is then clearly 
in the asymptotic regime as shown by comparing the curves 
with the symbols on the y-axis, which are the values obtained 
by directly solving the fixed-point equation. The two curves 
are superimposed, implying fjk ~ rfk, until the Larkin scale 
(marked by a circle), at which point they separate. However, 
the fixed point values are such that 2r; — f; is unambiguously 
strictly greater than zero. 



a whole domain of d and N, including the Ising version, 
for which instead fj — rj j^ and that near the boundary 
of this domain the exponents verify fj c^ rj ^ 2'q. 

We have pointed out that a potential caveat of the 
derivations of the two-exponent scaling, which imply that 
the random-field strength is not renormalized by fluctu- 
ations at criticality, is that they are blind to rare events, 
regions or samples that are known to play an important 
role in random-field systems. The arguments leading to 
two-exponent scaling however apply to random-field sys- 
tems with long-range correlated disorder. In such sys- 
tems, the bare variance of the random field is no longer 
given by Eq. ([2]) but is characterized by an exponent p: 



The long-range piece of the second cumulant of the ran- 
dom field is then not renormalized by fluctuations.—^— 
In this case, one therefore rigorously has that 2?] — fj = 
p. This is for instance correctly given by the heuristic 
derivation leading to Eq. ^ as the mean random field 
in a finite system scales as L'('^^py^. Similarly, a gener- 
alization of Eq. ((HI is now also valid. 
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FIG. 3: Ratio (2-q ~ fj)/r] as a function of d for the RFIM in 
the vicinity of doR — 5.1 from the NP-FRG. Above doR, the 
ratio is exactly equal to 1 as the stable fixed point is associated 
with the d — >■ d — 2 dimensional reduction. For d < doR, the 
stable fixed point is characterized by the existence of a cusp 
in the cumulants of the renormalized random field. The ratio 
then decreases continuously from 1 as one decreases d. As in 
Fig. 1, there is a finite domain of d for which 2r; — 77 7^ with 
no ambiguity. 



ht^{x)h''{y) 



Sf_cu 



(12) 



In the case of the short-range RFIM, the two-exponent 
scaling description is exact close to the lower and (more 
trivially) upper critical dimensions. It moreover appears 
to provide an extremely good approximation in between, 
with 277 — 77'^10~^or less. In the light of our results, this 
implies that computer simulations, even elaborate and 
extensive ones as the recent study in [ISj, will likely be 
always inconclusive due to unavoidable residual errors 4^ 
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